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ABSTRACT 

Based  upon  the  existence  and  uniqueness  of  a  solution  to  the  linearised 
Lundquist  equations  established  previously «  the  modified  energy  principle  for 
the  (^stability  of  a  confined  toroidal  plasma  is  rigorously  justified.  A 
variational  principle  is  developed  to  find  the  infimum  oF~^,  and  an  esti* 
mate  for  the  maximum  growth  rate  is  obtained.  The  results  are  also  extended 
to  a  diffuse  pinch  and  a  multiple  tori  plasma. 
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SIGNIFICANCE  AND  EXPLANATION 


In  this  report,  we  shall  justify  rigorously  the  so-called  modified  energy 
principle  for  the  d-stability  of  a  confined  toroidal  plasma.  Intuitively 
speaking,  a  plasma  is  called  o-stable  if  it  does  not  grow  faster  than 
exp(ot),  where  o  ^  1/t  and  T  is  some  characteristic  time  needed  for 
fusion.  The  modified  energy  principle  claims  that  the  plasma  is  o-stable  if 
some  energy  functional  is  nonnegative  and  unstable  if  otherwise.  We  develop  a 
method  to  get  an  upper  bound  for  o  and  the  maximum  growth  rate  for  the 
plasma  is  also  obtained.  Furthermore  the  results  are  extended  to  the  case 
that  the  plasma  fills  up  the  whole  conducting  shell  and  that  there  are  several 
plasma  tori  in  the  conducting  shell. 
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STABILITY  THEORY  OF  A  CONFINED  TOROIDAL  PLASMA 
PART  II*  MODIFIED  ENERGY  PRINCIPLE  AND  GROWTH  RATE 

Peter  I«aurence*  and  N«  C*  Shen^ 

§ 1 •  Introduction. 

The  MHD  stability  of  a  plasma  equilibrium  confined  in  a  magnetic  field  is 
one  of  the  most  iu^rtant  topics  in  controlled  thermonuclear  research.  The 
approach  to  this  problem  is  very  often  based  upon  the  linear  energy  principle 
formulated  by  Berstein  et.  al.  (1957).  In  the  justification  of  the  necessary 
condition  for  stability,  they  assumed  that  the  eigenfunctions  of  a  certain 
linear  operator  form  a  complete  orthonormal  basis.  This  assumption  may  limit 
the  scope  of  the  linear  energy  principle  in  applications.  Laval  et.  al.  (1965) 
relaxed  this  restriction  and  established  a  modified  energy  principle  for  the 
so-called  o-stabillty  of  a  confined  plasma,  including  the  linear  energy 
principle  as  a  special  case.  However,  in  their  derivation  the  existence  of  a 
classical  solution  to  the  corresponding  system  of  the  linearized  Lundquist 
equations  is  tacitly  assumed,  and  at  present  no  such  a  solution  is  knom  to 
exist.  In  this  report  we  shall  state  and  prove  rigorously  a  precise  version 
of  the  modified  energy  principle,  for  a  confined  toroidal  plasma,  which  is  the 
main  contribution  of  our  work.  Our  approach  relies  upon  the  existence  of  a 
generalized  solution  to  the  linearized  Lundquist  equations  we  established  in  a 
previous  report,  which  we  shall  refer  to  as  Part  I. 
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The  con ten te  of  this  report  are  briefly  explained  as  follows*  Zn  Section 
2,  we  present  the  rigorous  proof  of  a  precise  version  of  the  modified  energy 
principle  for  o-stability*  A  necessary  and  sufficient  condition  Is  given  for 
an  equilibrium  to  be  o-stable*  The  proof  of  the  sufficiency  essentially  fol¬ 
lows  the  results  obtained  in  Part  1}  that  of  the  necessity  needs  some  regular¬ 
ity  and  embedding  theorems  by  Lions  and  Nagenes  (1972 )•  In  Section  3*  a  vari¬ 
ational  principle  is  developed  to  find  the  infimum  of  o  for  o-stability. 

Zn  Section  4,  an  estimate  of  the  maximum  growth  rate  is  obtained*  Finally,  we 
extend  previous  results  to  the  cases  of  a  diffusim  pinch  and  a  multiple  tori 
configuration  in  Section  5« 
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|2.  Necessity  and  Sufficiency  of  the  Modified  Hhergy  Principle* 

We  begin  with  the  definition  of  o-stability  of  a  plasna  equilibriun* 

tefinitJ^qn  An  equilibrium  is  called  a-stable  if  all  solutions  (C(t),A(t)) 

of  the  evolutionary  variational  problem  EVP  (§6,  Part  I)  have  the  property 

that  there  exist  some  constant  c  >  0  and  o  such  that 

2 

■  Cl^  <  c  exp  2at  on  tO,**)* 

ZfV  — 

We  now  prove 
^eorem^  2*  1  * 

An  equilibrium  subject  to  the  conditions  of  |3,  Part  I  is  o*stable  if 
and  only  if 

a((e,A),(e,A))  +  ^  >  0,  V(5,A)  e  W.  (2.1) 

Proof; 

To  prove  the  sufficiency^  assume  that  A(t)  »  (C(t),A(t))  is  a  solution 
to  the  EVP  with  (€(0),A(0))  »  “  ^0  prescribed  flux 

/  A*nd8«o(t)  (|6,  Part  I)*  As  in  Part  1,  we  associate  with  the 
EVP  an  auxiliary  EVP  (AEVP)  for  which  we  construct  the  m'^"  Galerkin 
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We  let 

b  -  »«(0)l2,p  Ci*A(0)lJ  , 

and  would  like  to  pass  to  the  limit  in  (2.3).  Since  the  equilibrium  satisfies 

the  assumptions  in  §3,  Part  1,  and  by  the  coerciveness  of  a(*«*),  there 

2 

exist  A,  6  >  0,  (where  X  >  o  )  such  that 

a((5,A),(5,A))  +  Xiei!  ^  >  6|(C,a)|5  . 

Thus  proceeding  as  in  §6p  Part  I  for  the  existence  proof  and  using  the  renark 

in  §7p  Part  I,  we  have 

e“(t)  ♦  C(t)  in  L^[[0,T]  j  1.^(0  )]» 

P  P 

l“(t)  ♦  Ut)  in  L^[[0,T1  .  lJ(Q^)].  (2.4) 

P  P 

Also  given  the  construction  of  our  Galerkin  approximations#  we  have 

c“(0)  -*•  C(0)  in  L^{Q  ), 

P  P 

i“(0)  ♦  1(0)  in  1.5(0  ),  (2.5) 

P  P 

A“(0)  ♦  A(0)  in  W. 

For  TpS  e  [OpT] p  T  >  Sp  we  now  integrate  both  sides  of  (2e3)  from  a  to 


T  and  imiltiply  both  sides  by  1/(t  *  s)p  to  obtain 

8  8  8 

We  now  can  pass  to  the  limit  on  the  RelleS*  becuase  of  (2e4)  and  on  the 

LeHbSe  %«e  use  lul  <  lim  influ^l,  valid  for  any  weakly  convergent  sequence  in 
a  Banach  space#  to  obtain 

dt  «  J  .  f).^«dt  -  /  5  dt.  ,2.6, 

8  8  8 

Noting  that  ■C(t)l2  ^  6  I.^[0,T]  L^tO.Tj,  ws  now  1st  T  ♦  s  on  ths 
L.H.S.,  apply  ths  Lsbssqus  diffsrsntiation  thsorsm  [Rudin,  1974]  thsrs,  which 
statss  that  ths  dsrivativs  of  ths  intsgral  of  an  function  is  ths 

function,  and  by  explicit  integration  and  diffsrsntiation  on  ths  R.H.S.,  ws 


'4' 


obtain 


„  <  (15(0)1=  „ 

<  (15(0)1=  „  ♦  £).=”'. 

In  addition.  If  we  Insert  (2«7)  Into  (2*2),  we  have 

IC”(t)l2  p  <  +  X{(ie{0)l2  p  +  ^  > 

-  b"  +  ^(K(0)l2  p  -  b“. 

Upon  going  through  the  same  procedure  as  between  (2.5)  and  (2*6),  we  may  pass 
to  the  limit  and  obtain 

IC(t)l^  „  <  b  +  X(IC{0)l5  „  +  Y)e^®^  -  b 

-  X(I5(0)I^  p  +  (2.8) 

This  establishes  the  sufficiency  of  (2«1)« 

Next  we  prove  the  necessity  of  the  modified  energy  principle  for 

a--8tablllty»  We  saelc  to  show  here  that 

a((U),(i.A))  <  0 

P 

for  an  element  (C.A)  6  W,  implies  the  existence  of  a  solution  of  the  EVP, 
(C(t),A(t)),  such  that 

1^1 2  p  ^  c  exp  2at. 

In  the  proof  of  this  result  the  essential  part  is  played  by  certain  regu¬ 
larity  and  embedding  theorems  of  Lions  and  Nagenes  (1972),  which  are  then 
coupled  with  the  original  proof  of  Laval  et.  al.  (1965)  to  yield  the  results* 
He  begin  by  slightly  changing  our  point  of  view  and  noting  that  l(C,A)l^ 

defines  a  norm  on  P^  (H)  which  %re  will  denote  by  •  This  is  clear 

P  ^ 

since,  given  C  8  Q  ,  the  boundary  value  problem 
P 
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V  X  V  X  A  «  0 


V  •  A  •  0 


n  X  A  »  -(n  •  C)B^  on  T  , 

P 

n  X  A  -  0  on  F  , 

V 

/  A  *  n  ds  *  Of  (2.9) 

r 

V 

determines  A  uniquely,  and  the  dependence  on  C  of  A  is  linear. 


Let  Be  the  dual  space  of  W^.  He  seek  to  show  that  the  following 


formula  holds. 


2  2 
<  ^  e(t),e(t)>  ^  ^  <5(t),5(t)> 


3t' 


3t 


2,P 


(2.10) 


2  <3^  5(t),  3^  ' 


O  w 

where  <*,*>  ^  denotes  the  action  of  the  linear  functional  C(t)  e  H^ 

u  vu  at^  ^ 


WgXw^ 


on  C(t)  e  H^f  and  where  (C(t)fA(t))  is  a  solution  to  the  ASVP.  We  shall 

3 

also  show  that  <C(t)f$(t)>2  p  is  absolutely  continuous  on  [0,T].  How 
a(  (Cfk)  f  (Cf A) )  defines  a  bilinear  form  on  W^,  «diich  we  shall  denote  by 
a^(CfC). 

In  order  to  prove  (2.9),  note  that  since  (C(t),A(t))  is  a  solution  of 


AEVP 

32  ^ 

<  —Z  C(t),?(t)>  ^  +  a.(?,5)  -  0,  v4  e  w  -.  (2.11) 

at^  w^xw^  ^  ^ 

We  now  invoke  a  theorem  of  Lions  and  Nagenes  (1972).  The  solution  of 

(2.11)  has  the  property 

5(t)  e  c°{ [o,T]  t  w^), 

a£  0 

*5“  e  c  (to,T]  I  Hg)  , 

2 

where  denotes  the  closure  in  (^p)  of  W^,  previously  denoted 

Pq  (H^).  Furthermore,  there  exist  C^(t)  e  C^((0,T]  :  W^),  such  that 

a  e  1 

e  (t)  e  C  ([0,T]  i  W^)  and 
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(2.12) 


ret)  ♦  c(t) 


in  c"([0,Tl  S  w^), 


It  ■*’11^*'^  C°([0,T]  :  Hg)  , 

3^5®  3^5  2  * 

(t)  + -“(t)  in  L‘‘(tO,Tl  :  W.)  . 

3t^  3t^  ^ 

We  shall  first  assume  that  formula  (2.10)  holds  for  C 

holds  for  C.  Integrating  (2.10)  from  0  to  T,  t  6  (0,T] 


(2.13) 

(2.14) 
euid  show  that  it 

we  have 


r  2<  ~  c®(t),c®(t)>  . 


at" 


W^xw^ 


-  2  /^<  C^(t),  •It’  5^(t)>,  _dt. 


(2.15) 


■St  ^  '■"2,P 

Now  using  (2.13)  we  may  pass  to  the  limit  an  the  integral  on  the  L.HmSm  of 
(2m14),  and  using  (2*11)  and  (2*12)  we  may  pass  to  the  limit  in  the  terms  of 


the  R.HwSm  Thus  we  obtain 

.2 


r  2<  C(t),e(t)>  ^ 


at' 


Tt 


<C(t),C(t)> 


t»T 
2,p|t-0 


-  2  /^  <  Ht),  C(t)>2  pdt. 

Again  since  all  integrands,  because  of  (2.12)  to  (2.14),  are  in  L^[0,T] 

thus  in  L^(0,T],  Lebesgue's  differentiation  theorem  allows  us  to  differentiate 

the  integrals  on  both  sides  of  (2»16)  with  respect  to  t  and  recover  the 

a 

integrands,  thus  it  follows  that  (2.10)  holds  and  that  <C(t),C(t)>2  ^  is 
an  absolutely  continuous  function  on  [0,T]  being  equal  to  the  integral  of 
an  function. 

We  now  return  to  show  that  (2.10)  holds  with  $^(t)  replacing  C(t). 

For  this  purpose  we  use  an  embedding  theorem  of  Lions  and  Magenes  (1972)  which 
reduces  to  the  follo%ring  special  fora  needed  later. 
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? 


/ 

i 


r 


If  n(t) 

has  the  properties: 

1) 

n(t)  8  L^({0,TJ  ;  W^) 

2) 

e  I.^((0,TJ  :  W^), 

3^1)  2  ,  ,  * 

3) 

- ^  (3  L  (I0,T1  :  W-) 

3t  ^ 

then  n(t)  6  A  C( [0,T]  i  W^), 
e  AC([0,T]  ;  H^). 

Moreover,  If  we  let 

X  =  {ulu  6  L^([0,T]  ;  W^),  ||-C  L^(10,T]  ;  W^), 

3  2  * 

— ^  e  L^([0,TJ  ;  W  )  , 

3t^  ^ 

with  norm 


lul  -  {lul%  ,  + 

,9^u.2 


W^) 


3t^  L^([0,T]:W^) 


(2.18) 


(2.19) 


then 

Cq((0,T]  :  W^)  is  dense  in  X* 

m 

Thus  in  particular  C  ([OrT]  :  W^)  is  dense  in  X« 

g 

Given  that  C  (t)  satisfies  the  above  conditions  (cf«  (2.12)  to  (2*14)), 
vie  nay  take  (C^)^(t)  6  (C  ((0,T)  :  W^)  >  5^<t)  in  X.  Now  starting  with  the 
relation 


3^(5®) 


r  2<  - 5-2- (t),  (?®)„(t)>  ^  dt  -  ■|-<(5®)„(t),(5®)„(t)>- 

0  2t2  "  w*x  Wg 


t«T 
P  t«0 


-  2  r  <  It  It 


for  T  6  [0,T],  and  using  (2.17),  we  may  pass  to  the  limit  in  the  above 


equation  with  the  result 
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2  C 

r  2<  (t),  ^  ^  <c®(t),5®u)> 


at" 


2,P 


t*T 

t*0 


-  2  /^<  C*^(t), 


for  t  9  [0,T],  So  we  obtain  (2,15)  as  desired.  This  concludes  the  proof 
that  (2.10)  holds  with  C^{t)  replacing  C(t).  We  now  invoke  another  lemma 


of  Lions  and  Magenes  (1972)  which  establishes  that  the  energy  equality 
3£  2 

holds  for  solutions  (C,A)(t)  to  AEVP. 

The  proof  now  proceeds  almost  identically  to  that  of  Laval  et.  al« 

(1965).  For  completeness  sake  we  present  it  here. 

If  C(t)  is  a  solution  of  the  AEVP  defined  by  (2*11),  then 

<  — r  C(t),e(t)>  .  +  a.(C,5)  =  0 

at^  w^xw^  ^ 


(2.20) 


(2.21 ) 


holds  because  C(t)  €  W^. 

If  we  make  use  of  (2.9)  and  substitute  for  a(C,C)  from  (2.19)  into 


(2.21),  we  obtain  the  virial  equation 

We  specify  the  initial  data  for  AEVP  to  be 


(  I4i 


o?  . 


(2.22) 


(2.23) 


1 1*0  "  3 1 1 1=0 

Note  that,  with  this  choice  of  initial  data,  C  in  (2.22)  is  negative,  so  %ie 


obtain 


on 


[0,T). 


Thus  the  same  argument  as  from  (6.20)  to  (6.21),  Part  I,  shows  that 

n  on  [0,T] 

2/P  ^^2  2,p  dt  2,P 


(2.24) 

Note  also  that,  since  according  to  (2.19)  C(t)  9  AC((0,T]  :  H^),  and 
since  IC(0)I  «  Kl  ^  0,  if  we  choose  T*  <  T  small  enough.  >  0  such 

that 

IC(t)l.  ^  >  a  for  t  9  [0,T*).  (2.25) 
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Thus  in  particular  if  we  define 

•  r  #  ^  2 


ie(ty-: 


y(t)  -  log 


y(t)  is  absolutely  continuous,  since  the  logarithm  is  a  Lipshitz  function  for 
any  possitive  interval  bounded  away  from  zero,  and  since  the  composition  of  a 
Lipshitz  function  with  an  absolutely  continuous  function  is  absolutely  oontin-* 
uous  ( S  tampachi a ,  1 9 66 ) • 

Similarly, 

is,  because  of  (2.19),  the  quotient  of  an  absolutely  continuous  function  and 
an  absolutely  continuous  function  which  is  bounded  away  from  zero,  and  hence 
(Rudin,  1974}  is  absolutely  continuous* 

Now  (2*24)  becomes 

>0  on  [0,T' ), 

dt^ 

with  the  initial  conditions  y  *  0,  ^  *  20.  Thus,  since  y(t)  and 


^6  AclO,TM,  we  may  integrate  this  inequality  and  obtain 


^  >  20,  and  y(t)  >  20t  for  t  G  [0,T"], 


(2.25) 


therefore. 


I5(t)l^  n  >  •5(0)1^  „  exp  20t  for  t  6  t0,T']. 

^  fl'  *  ,P 


(2.26) 


We  now  note  that  the  condition  (2.25)  is  clearly  satisfied  for  all  t  >  0 


by  repeating  the  same  argument,  we  obtain 

n(t)il  „  >  IC(0)|J  „  exp  20t. 
2,p  2,P  ^ 

This  concludes  our  proof  of  Theorem  2.1. 


S3*  Variational  Principle 


Earlier  we  established  that  a( (C,A), (C,A))  is  coercive  with  respect  to 
the  W  norm  on  W  x  w*  There  exist  31  >  0  and  6  >  0  not  necessarily 
unique  such  that 

a((e,A),(e,A))  +  Aicij  >  6I(C,A)I^,  V(C,A)  6  W. 

In  particular,  X  may  be  chosen  so  that 

a((C,A),(C,A))  +  p  >  0  (3.1) 

V(e,A)  6  W.  (3.2) 

We  now  show  how  to  find  the  minimal  X,  X^,  that  ensures  (3*2)  by  a 
variational  procedure* 

We  wish  to  find 


inf 

ntl  =1 

2,P 


a((e,A),(C,A)). 


(3.3) 


(C,A)ew 

As  seen  from  (2.1),  this  Infimum  Is  finite.  We  now  show  that,  although  in 

general  there  is  no  element  (C,A)  6  W  which  realizes  the  infimum,  we  can 

construct  a  sequence  of  vectors  (C”»a")  e  W  such  that 

a((t'‘,A",(C",A"))  -  X  +  X^  . 

n  u 

We  use  again  the  skew  orthonormal  basis 


<  >2,P 


'ij' 


^  n  ^  ^ 

and  define  A  »  (€  ,A  )  to  be  an  element  of  W  of  the  form 
4”  .  (?",»»,  . 


t3.4) 


with  1C  •_  „  =  1,  for  which  the  infimum  of  (3.3)  is  attained.  That  this 
P 

infimum  is  indeed  reaalized  by  an  element  of  w  of  the  form  (3*4)  follows  by 
noting  that 


■  il  • 
■5"'"  -  I  •?.  . 


80  that  the  vary  over  a  conpact  set  of  dieension  n«  Since 

a(  is  a  continuous  function  of  a  ,  it  attains  both  its 

in 


maximujB  and  minimiiai  and  the  ninimua  satisfies 
a(A",P^)  -  1  <  i  <  n, 

a"  .  a(A",A"). 

The  equation  (3.5)  holds  because  Vr  >  0,  v£  of  the  fora  (3.4), 


(3.5) 


(3.6) 


£  «  such  that 

a(A*‘  +  rE,A"  +  rl)  >  +  ri:,,C"  +  E.>,  ^ 

1  1  2,P 

a(A",A")  +  r^a(E,E)  +  2ra(A",E)  >  x"<i",cV  „  +  xV<E,,r  >  ^ 

2fP  1  2  29P 

+  2X"r<Z^,C>2  p. 


Using 


a(A",A'')  -  X"<C",C">  -  X"  , 


to  cancel  one  term  from  both  sides  and  then  dividing  by  r  and  letting 
r  ♦  0,  we  have 

a(A",r)  < 


Therefore 


a(A",E)  >  ^"<5",E,>2^P 


follows  from  replacing  £  by  >£  In  the  above  derivation.  Now  X**  Is  an 
Increasing  sequence  of  real  numbers  since  the  Inflmua  Is  ta)cen  over  an 
Increasing  sequence  of  sets.  Moreover,  by  the  coerciveness  there  exist  X 


such  that 


a(A'',A'')  +  Xie"!,  „  >  0, 
2#P 


that  is. 


X  >  -X. 
n 


Hence  the  sequence  of  {X^}  is  bounded  below  and  there  exists  X^  such  that 


where  X^  may  be  negative.  Moreover,  by  the  coerciveness  again. 


6,a  >  0 


such  that  a  > 


a((5,A),(C,A))  +  >  61(5, A)l^  , 


(3.5) 


where  n  "  follows  that 

^ ,  P 

'In  "n 
A  =  ,A  ), 


forms  a  bounded  sequence  in  W,  since 


a((5",A"),(5",A'‘))  +  >  Xq  +  a"  >  6l  (5'',a'')  1^^, 


Vn. 


Hence  there  exists  a  weakly  convergent  subsequence  in  w  of  the  A  ,  and  in 


particular,  of  the  in  L^(fl  )< 

P  p 


Unfortunately,  without  some  additional 


features  of  w  ensuring  compactness  of  the  embedding  of  W  in  L  (ft^)  x 

2  ^ 
h  we  may  not  conclude  that  the  limit  A  is  a  genuine  eigenfunction. 

We  will  not  examine  this  question  further  here.  In  any  case  we  have  shovm 

a(A",A")  =  X"  +  +X„ 


where  a”  is  the  solution  to  the 


minimization  problem  discussed 


above.  Finally  we  have 
Theorem  3. 1 . 

inf  a( (5, A) , (5, A) )  »  X 

ml  =1 

2,P 


Proof; 


(5,A)6W 


inf  a( (5,A) , (5,A) )  <  X  , 

ml  =1 

2,P 

(5,A)ew 

A  A 

is  clear  since  (4  nA  )  is  a  collection  of  elements  of  W  with 


'2,P  = 
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inf  >  X 

•'•Ip- 


(C,A)6W 

is  shown  by  contradiction p 


If 

inf  a{(C,A),(C,A))  <  X  -  e 

■'^•Ip- 

(C,A}eif 

then  there  exists  H  >  0  such  that  for  n  >  N 

ad",!**)  <  Xjj  -  ^  . 

But  this  contradicts  the  defining  property  of  A"  and  the  Mnotonicity  of  the 
sequence 


A  consequence  of  (2«7)  is  the  followingt 
Corollary; 

a((C,A),(5,A))  -  ^  >  0 

0  2,p 

V(€rA)  6  If. 


Proof: 


Since 

inf  a((5,A),(5,A))  -  X  , 

2»P 

(CpA)a« 

a((C,A),(C,A))  -  Xq  >  0. 

V(C,A)  6  W  with  151?  „  -  1. 

2pP 

Now  for  arbitrary  (5/A)  e  N,  consider 

(5,A)  -  (5,A)/I5I^  -  , 

where  ^  •  Then 

P  P 

a({5,A),(5,A))  -  Xq  >  0. 

However,  due  to  the  bilinearity  of  a(*,*) 

a((5,A),(5,A))  -  a{(5,A),(5,A))/l5lj  p. 

Ne  obtain  (2.9)  by  multiplying  (2.10)  by  15I2  p* 
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§4.  Estimate  of  the  Maximum  Growth  Rate. 

In  deriving  the  sufficiency  of  the  modified  energy  principle  in  Theorem 

2 

2.1,  we  established  that  if  ao  such  that 

a((C,A),(C,A))  +  >  0,  V(C,A)  e  w, 

then  any  solution  (C(t),A(t))  of  the  EVP  on  (0,*)  is  subject  to  the 
growth  estimate 


p  —  ®  to,**).  (4.1) 

In  the  process  of  proving  the  coerciveness  of  the  bilinear  form 
a( (C/A) (Cf A) )  (fS,  Part  I),  that  is,  aX,6  >  0  such  that 

a((C,A),(C,A))  +  ^»e*2,P  ^  va,k)  6  H 

we  found  that  upper  bounds  on  A  could  be  given  in  terms  of  the  constants 
e^,  which  characterize  the  equilibrium  quantities  p,  B,  V  x  b.  However,  in 
|5,  Part  I,  6  on  the  R.H.  S.  of  equation  (5.16)  was  taken  to  be  V2  •  now 
make  note  of  what  estimate  is  obtained  for  A  when  6  »  0,  since  this  is  all 
that  is  needed  in  order  to  obtain  the  growth  estimate  (4.1  )• 

Thus  tie  let  f^  and  £2  in  (5.15),  Part  I  be  such  that 

«l£l  -  ^  ®3®4^2  “  ’ 


X 


•l*2 


+ 


4 


is  an  upper  bound  for  the  naximui  growth  rate;  or  we  aake  use  of  the  w>re 


precise  estiaate  following  froa  (5>13)  and  (5.14)  in  PartZ, 
,  .2V  X  B  >«  n  •  B  •  ViK 

•  » 


is,  again,  an  upper  bound  on  the  aaxiaua  growth  rate. 
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§5*  Diffuse  Pinch  and  Multi -Tori  Configuration* 

In  this  section  we  extend  the  results  of  $3,  Part  I  to  two  other 
configurations  of  importance  in  thermonuclear  research*  One  case  deals  with  a 
diffuse  pinch  and  the  other,  a  multiple  tori  plasma  confined  in  a  toroidal 
shell* 

In  a  diffuse  pinch  (Figure  1)  the  plasma  extends  to  the  confining 
shell*  The  equations  in  the  plasma  region  are  the  same  as  in  ^3,  Part  I;  now 
we  only  need  to  express  them  in  Eulerian  coordinates 

as 

PC  -  V(YPV  •  5  +  e  •  Vp)  +  V  X  B  *  V  X  (C  *  B))  (5.1) 

+  7  X  y  X  (g  X  B)  X  B, 

We  impose  the  boundary  condition  on  the  confining  shell  that  the  normal 
component  of  velocity  on  the  shell  is  zero.  By  integration 

^•n>Oonr,  (5.2) 

and  we  note  that  if  (5.2)  is  satisfied  then  the  perturbed  magnetic  field 
B^  *  V  X  (g  X  B) 

automatically  satisfies  the  boundary  condition: 

B^  •  n  ■  0  on  P, 

because 

g  •  n  -  0  »>  n  x(g  X  B)  -  0  on  P, 
where  B  *  n  ■  0  by  our  assumption  on  the  equilibrium  state/  and 
n  X  V  >  0  •■>  n  •  7  X  V  ■  0  on  P, 

as  CM  be  easily  seen  by  applying  Stoke 's  theorem  to  infinitesimal  loops  on 
the  boundary. 

We  now  present,  in  a  somewhat  condensed  fashion,  m  existence  and 
uniqueness  proof  for  a  weak  solution  to  the  boundary  value  problem  (5.1)  to 
(5.2)  with  initial  conditions 
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(5.3) 


N«  shall  Mke  the  following  assuaptions  on  the  equilibrlua  quantities t 


P,  Vp  e  L  , 

B  ■  n  «  0  on  r. 

Furtheraore,  we  also  assuae  ac^,C2*d^,d2  >  0  such  that 
<  p(x)  <  Cj  , 

d,  <  P(x)  <  d^. 

If  we  take  the  inner  product  of  (5.1)  with  C  where  K  satisfies  the 

boundary  condition  (*n>0  on  T,  we  get 

«  ■  *  /  •  CV  •  C  +  V  X  (5  X  B)  •  V  X  (C  X  B) 

2,9  Q 

P 

-5»Vxbx  (7x  (^xb))-!-  (V»C)(5*  Vp)dv 

I  /  (C  •  n)lTPVe  -B»  {7x5xB)-e*  VPldS 
'p 

-  -  /  YP(7  •  0(7  •  ^)  +  7  X  (5  X  B)  •  7  X  X  B) 

-  6  *  7  X  B  X  (7  X  (C  X  B))  +  (7  •  0(C  •  7p)dv. 

Thus  ws  define 

a((C,?)  -  /  YP(7  •  €){7  •  O  +  7  X  (5  X  B)  •  7  X  (5  X  B) 

“P 

-  5  •  7  X  B  X  (7  X  (5  X  B))  +  (7  •  0(C  •  7p)dv,  (5.4) 

and  obtain 

<pi>C>-  _  +  a(C>0  *0  vC  such  that  5  •  n  -  0  on  T, 

2pP 

Let 

v"  -  (C  e  H^(0)|(C  •  n)  -  o}, 
with  scalar  product 

<5,0  -  /  (Y.(7  •  5)(7  •  5)  (5.5) 

“P 

+  7  X  (C  >c  B)  •  7  X  (5  «  B)  ♦  p5  •  ?)dVp 
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and  let 


V  »  closure  of  V  in  I  I  • 

v" 

Ml  the  properties  related  to  volume  integrals  enjoyed  by  the  first 
components  of  elements  of  W  (cf.  §4,  Part  I)  also  hold  for  elements  of  V 
and  the  proofs  are  the  samef  that  is 
7  •  c  e 

7  X  (C  X  B)  e  L^(a), 

5  e  L^{Q). 


We  now  make  use  of  a  result  in  Temam  (1973),  noting  that  given  (5«5), 
3c  >  0  such  that. 


so  that  if  is  a  Cauchy  sequence  in  V*, 


♦  C  in  V, 


with 


then 

n  •  ♦  n  •  5  in  h"  ^2  (d. 

Therefore,  since  n  •  €*  *  0  in  H*  ^2  (P)  (in  fact  in  H  ^2(r))), 
n  •  C  »  0  in  h"  (T). 

All  Other  aspects  of  the  existence  and  uniqueness  proofs  are  straight¬ 
forward  simplifications  of  those  in  |6,  Part  X,  so  we  will  not  repeat  them 
here  but  only  state  the  follo%#ing 


Theorem  S«1; 

There  exist  a  unique  solution  C(t)  of  the  EVP 

(1)  ^  <C(t),C>.  ^  +  a{C,?)  -  0  VC  e  V, 
a^t 

(2)  5(0)  -  5q  e  V,  ||(0)  -  e  l^{0), 
and 

5(t)  e  C®(10,«)  :  V), 

11  e  C°{(0,«)  :  L^(0)), 
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One  cen  easily  verify  that  for  the  diffuse  pinch  the  estleate 

Q 

IC(t)lJ  „  <  tn  +-^5(0)1^  ^  +  ^C(0)l^)e^®^ 

2,0  g2  2,p 

holds.  Where 


2 

o  can  be  estimated  from  (5.13),  Part  I, 

2pxBXn*B*V 
-2  .  n. 

o  .  I - - - 1 

L 

If  the  plasma  region.  Instead  of  consisting  of  one  torus,  is  the  union 
of  N  non-intersecting  tori  with  interior  Qp^,  all  enclosed  in  an  outer 
toroidal  conducting  shell  (Figure  2)  F^,  the  basic  approach  of  Part  1  again 

applies. 

In  each  connected  plasma  domain  Qp^  we  search  for  C^(t)  such  that 

P^C^(t)  -  P(5^)  -  V(YpV  •  5^  +  •  Vp^)  (5.6) 

+  7  X  X  V  X  (C^  X  B^)  +  7  X  7  X  x  B^)  x  b^. 

In  the  vacuun  region  we  gain  require 
V  X  V  X  A  «  0, 


V  X  A  »  0. 

In  each  plasna-^vacuua  interface 


the  continuity  of  the  total 


pressure  requires  that 

,  .  t  .  i 

-  7p  V  •  +  B^  (5  •  ^  +7 

-  B^  •  (7  X  A  +  5^  •  VB  ''), 

PurtlMrmore, 

n  X  A  -  0  on  .  (5.10) 

We  also  have  the  prescribed  flux 

/  A  •  n  ds  •  6(t),  (5.11) 

r 

V 

and  initial  conditions 


X  (r  X  B**  )) 
on  r  ,  . 


(5.9) 


C(0) 
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Proceeding  in  the  sane  way  as  in  Part  I,  we  define  a  space 


w”  »{(€’,•••, e  h’(Qp^),  a  e  v^CQ^)  h’(Q^)} 


n  X  A  >  0  on  r  , 

V  ' 

n  X  A  ■  (-n  •  on  F 


with  l(CrA)l  defined  by 

w“ 


i.2 


(5\5^,**-,C**,A)I  -/  |v  X  A|^dv  +  X  /  {ifP^{V  •  ch 


H  Q 


i  np4 


+  |V  X  (€^  X  B^)|^  +  p^[}dv. 


and  let 


W  *  closure  of  w"  in  I  I 


W 

w'*’  -  W  •  (0,A'),  A*  8  Vj  , 

We  note  that  in  order  to  determine  a  unique  solution  to  the  elliptic  system: 

V  X  V  X  A  -  0, 

V  •  A  »  0, 

n-1 

n  X  A  -  0,  on  [  U  F  ]ur  , 

1-1  Pi  ^ 

it  Is  now  necessary  to  prescribe  [Blank,  Grad,  Predrlchs,  1957]  N-1  fluxes 

/  A  •  n  ds  «  6. (t)  , 

'■'>1 

on  any  M  1  of  the  N  plasma  surfaces  say  F  ,***pF  ,  in  addition  to 

Pi  Ph-1 

(5.11 ). 

With  conditions  on  the  equilibrium  quantities  identical  to  those  of  |3, 
Part  Ip  and  the  same  procedures  as  in  §3  to  §6p  Part  Ip  we  may  again  establish 
the  existence  and  uniqueness  of  solutions  to  the  following  ABVP  associated 

with  (5.6)*(5.11): 


22* 


(1)  “  <5, (5, - 

3t^  1  1  z,p  ^ 


(5,12) 


v((r, A)  e  w. 


(2)  (5\o),***,c”(0),A(0))  - 

(||-(0),‘*«,||-(0))  =  e  H, 


(5.13) 


where  H  »  (closure  in  L  of  ^  solution  of  (6«13), 


corresponding  to  fluxes  6(0)  and  6^(0). 


(3) 


/  A  •  n  ds  *  0,  /  A  *  n  ds  4  0. 

r  r 


(5.t4) 


Pi 


Theorem  5»2; 

There  exists  a  unique  solution  of  the  AEVP  (5.1 2-5.1 4)  such  that 


(c’(t),'*‘,£**(t),A(t))  e  w  , 

(|^t),*«»,|^t))  ec°([0,«)  ;  H). 

To  the  solution  o£  AEVP  we  add  a  solution  of  the  elliptic  system 
7  X  7  X  A*  =  0, 

7  •  A'  -  0, 


n  X  A'  »  0  on  V  rp.  U  r  , 

i  1  V 

/  A*  •  n  ds  *  6(t),  /  A*  •  n  ds  »  ^j(t), 

r  r 

'L  '’i 

where  A*  8  V^. 

As  before  this  solution  has  the  property  that 
2 

(1)  “  <(5\***,e”)(t), (?’,•••, 5“)>,  „ 

at 

+  a((c\«»*,c“)(t),A(t)(«\*-*,C**(t),A(t)) 


v(l\***,l*'(t),A(t))  e 


where 


0 
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«  I  a{(r,A^),(r,A^)). 

i»1 

(2)  ((5\o),‘**,C*'<0)),A(0))  .  ((eo,*“,eJ).AQj  g 

1  N 

(§|-(0),*-,||-(0))  »  (CQ,***,eJ)  e  H. 

(3)  /  A  •  n  ds  “  ®(t);  /  A  •  n  ds  »  a. (t). 


rp 


If  o(t)  is  continuous,  then 


ij 


{(€’(t),--*,C**(t)),A(t))  e  C{(0,T]  ;  w"^) 

3£^  3£**  0 

e  C  {[0,T1  I  H). 


Here 

It  can  be  shown  that  growth  estimates  of  the  form 

15^1?  „  <  C  exp  2ot 
2 

again  hold,  tdiere  o  is  the  smallest  positive  real  number  such  that 


i-N 


a((5\***,e**,A),(c\**«,5",A))  +  I 


i-1 


2,P 


>  «l(c\«‘*,c"),A)l^. 
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tional  principle  is  developed  to  find  an  infimum  of  0,  and  an  estimate  for  the 
maximum  growth  rate  is  obtained*  The  results  are  also  extended  to  a  diffuse 
pinch  and  a  multiple  tori  plasma* 
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